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ABSTRACT 

Nuclear receptors (NRs) are central regulators of 
pathophysiological processes; however, how their 
responses intertwine is still not fully understood. 
The aim of this study was to determine whether 
and how steroid NRs can influence each other's 
activity under co-agonist treatment. We used a 
unique system consisting of a multicopy integration 
of an estrogen receptor responsive unit that allows 
direct visualization and quantification of estrogen 
receptor alpha (ERa) DNA binding, co-regulator re- 
cruitment and transcriptional readout. We find that 
ERoc DNA loading is required for other type I nuclear 
receptors to be co-recruited after dual agonist treat- 
ment. We focused on ERot/glucocorticoid receptor 
interplay and demonstrated that it requires steroid 
receptor coactivators (SRC-2, SRC-3) and the 
mediator component MED14. We then validated 
this cooperative interplay on endogenous target 
genes in breast cancer cells. Taken together, this 
work highlights another layer of mechanistic com- 
plexity through which NRs cross-talk with each 
other on chromatin under multiple hormonal stimuli. 

INTRODUCTION 

Estrogen receptor alpha (ERa) is a ligand-activated tran- 
scription factor that plays important pathophysiological 
roles in multiple tissues, including mammary gland, 
uterus (1,2) and bone (3,4). ERa is a member of the 
nuclear receptor (NR) superfamily characterized by 
modular structural domains (A-F) (5) that include trans- 
activation AF1 (A/B) and AF2 (E/F) regions, DNA- 
binding (C) and hormone-binding (E/F) functions. 
Ligand-bound ERa dimerizes and dynamically binds to 



estrogen response elements (EREs) throughout the 
genome; subsequently, ERa can serve as a platform for 
the recruitment of histone-modifying enzymes and tran- 
scriptional co-regulators that modulate chromatin struc- 
ture and RNA polymerase II action, impinging directly on 
target gene regulation. Co-regulator proteins, such as 
pl60s (SRC-1, -2 and -3) (6) and mediator (7) families, 
have been shown to greatly affect ERa-mediated tran- 
scription through direct interactions primarily with the 
AF2 region. 

Recent studies highlighted a varied and intricate web of 
different transcription factors that are responsible for 
finely tuned ERa-mediated transcriptional output. For 
example, ERa regulates the expression of other transcrip- 
tion factors, such as GATA3, FOXA1, MYC, PGR, FOS 
and PITX1 (8-11), which can then feedback on ERa sig- 
naling by working as pioneer factors, by marking the 
genomic landscape, by changing ERa expression and/or 
by directly modulating transcriptional output of second- 
ary target genes, thus creating a fine tuned phenotypic 
response to hormone. An additional layer of complexity 
resides in the cross-talk between ERa and other 
transcription factors when these are simultaneously 
activated by their own agonistic stimuli. This phenomenon, 
despite its physiological relevance, is less well studied. 
The best example of this in the NR field has been the 
analysis of the interplay between ERa and retinoic acid 
receptors in breast cancer cells, which showed both 
cooperative (12) and antagonistic (13) effects on ERa- 
mediated transcription after combination of estrogenic 
and retinoid stimuli. Other systems have shown similar 
interplay between NRs and transcription factors (14,15). 

Of the 48 NRs (16), seven comprise type I steroid 
receptors [e.g. androgen receptor (AR), ERa, ERp, pro- 
gesterone receptors (PR-A and PR-B), MR and gluco- 
corticoid receptor (GR)], and this subfamily is 
particularly important in a wide range of normal 
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physiological and pathological conditions (17,18). In the 
current study, we aimed to investigate the molecular 
cross-talk between type I NRs by using a unique tool de- 
veloped in our laboratory, the GFP-ERa:PRL-HeLa 
array cell line (19-22). The GFP-ERa:PRL-HeLa model 
consists of multicopy integration of an estrogen responsive 
transcriptional reporter gene derived from the regulatory 
region (promoter and enhancer) of the rat prolactin (PRL) 
gene [(23); e.g. PRL array] and also stably expresses 
GFP-ERa (green fluorescent protein-tagged) at levels 
comparable with MCF-7 breast cancer cells. The PRL 
array model allows for direct quantitative visualization 
of multiple transcriptional mechanisms via multiplex 
fluorescence microscopy, and it is readily amenable to 
the type of high-throughput analysis required to study 
complex systems level issues in gene regulation. 

After an initial immunofluorescence analysis of all 
type I NRs (minus MR), we decided to focus on the 
cross-talk between ERa and GR. It is increasingly clear 
that understanding the mechanisms underlying this 
cross-talk has great potential in studying multiple 
disease states, including inflammatory diseases and 
cancer. For example, in breast cancer, ERa, which is ex- 
pressed in ca. 70% of tumors (24), is currently the primary 
target for therapeutic intervention. Additionally, gluco- 
corticoids have been used for many years in the treatment 
of breast cancer to fight inflammatory processes and side 
effects of chemotherapy (25) with a wide range of effects 
observed in vivo (26,27) and in vitro. 

Enabled by quantitative analysis of GFP-ERa:PRL- 
HeLa cells using automated multiparametric image 
analysis, we show here that other type I NRs, but not 
other members of the superfamily, are recruited to the 
ERE-rich array only when ERa is present and under 
dual ligand treatment (e.g. in the presence of agonists 
for both receptors). Deletion mutant analyses reveal that 
co-recruitment of GR to the PRL array is dependent on 
the co-regulator interaction domain of ERa (e.g. helix- 12 
in the ligand-binding domain). Through a small co-regu- 
lator siRNA survey, we identified SRC-2, SRC-3 and 
MED 14 as essential ERa-binding partners required for 
GR recruitment. Importantly, fluorescent in situ hybrid- 
ization (FISH) studies indicate that co-recruitment of GR 
has a negative impact on ERa-mediated transcriptional 
stimulation. Next, we extended our analysis in breast 
cancer cells and found examples of genes and enhancers 
where the ERa/GR cross-talk occurs and has various 
gene-specific outcomes (e.g. cooperation and antagonism). 
Taken together, this study reveals unappreciated func- 
tional chromatin interplay between ERa and GR that 
depends on co-regulator recruitment and results in both 
cooperative and antagonistic transcriptional regulation as 
part of the fine tuning of hormonal responses. 

MATERIALS AND METHODS 

Cell culture and immunolabeling 

GFP-ERa:HeLa-PRL cells were maintained in Dulbecco's 
modified Eagle's medium containing 10% fetal bovine 
serum (FBS) (Gemini Bioproducts), 1 nM Tam (Sigma), 



200 ug/ml of hygromycin, 0.8 ug/ml of blastocidin, sodium 
pyruvate and L-glutamine. Cells were plated for 48 h in 
5% stripped-dialyzed FBS, phenol red free-Dulbecco's 
modified Eagle's medium before treatment on standard 
12-mm glass coverslips or 3 84- well optical glass 
(Greiner) or plastic (Aurora; siRNA experiments) plates. 

For immunolabeling, cells were fixed in 4% formalde- 
hyde in PEM buffer (80mM potassium PIPES, pH 6.8, 
5mM EGTA and 2mM MgCl 2 ), quenched with 0.1 M 
ammonium chloride and permeabilized with 0.5% Triton 
X-100. Cells were incubated at room temperature in blotto 
(5% milk in lx TBS-Tween 20) for 1 h, and then specific 
antibodies were added overnight at 4°C. The antibodies 
used were anti-GR (GeneTex GTX101120), RNA poly- 
merase II (ABCAM ab5408), PPARy (ABCAM 
ab52270), ERa (Millipore 04-820), Ser5-phospho RNA 
polymerase II (ABCAM ab5401), MED14 (Bethyl 
A301-044A), SRC-2 (BD Transduction Labs # 610985) 
and SRC-3 (BD Transduction Labs # 611105). After 
washes, the cells were incubated at room temperature 
for 1 h with Alexa Fluor 647 goat anti-rabbit IgG or 
Alexa Fluor 647 goat anti-mouse IgG (Molecular 
probes). This was followed by two PEM washes and 
another fixation and quench step before DAPI (1 ug/ml) 
staining. Coverslips were mounted on slides using 
SlowFade Gold (Invitrogen S36937). 

Fluorescence in situ hybridization 

Cells were fixed in 4% paraformaldehyde in RNase-free 
phosphate-buffered saline for 15 min and then perm- 
eabilized with 70% ethanol in RNase-free water at 4°C 
for 1 h. Cells were washed in 1ml of wash buffer (2x 
SSC (Ambion) and 10% formamide) followed by hybrid- 
ization in hybridization buffer (1 g of dextran sulfate, 1 ml 
of 20 x SSC buffer, 1ml of formamide and 8 ml of 
nuclease-free water) with RNA probes (dsRED2 
Stellaris™ probes, Biosearch Technologies Inc.) for 4 h 
at 37°C followed by one wash in wash buffer for 30 min at 
37°C and then DAPI stain for 10 min at 37°C. Finally, 
cells were washed in 2x SCC buffer followed by 
mounting using SlowFade Gold. 

Transfections, plasmids and RNAi experiments 

GFP-ERa, GFP-ERa (amino acids 1-534), GFP-ERa 
(amino acids 1-554) and GFP-GR plasmids were 
described previously (22). mCherry-GR was created by 
removing the GR from GFP-GR, using BamHI and 
BsrGI restriction sites and inserting it into the mCherry- 
C2 vector. Cells were reverse transfected using 
Lipofectamine 2000 reagent (Invitrogen 11668-500) for 
4-6 h. For RNAi experiments, we used Lipofectamine 
RNAiMax (Invitrogen 13778-150) along with 8nM 
siRNA. PRL-HeLa cells were transfected for 72 h. 

Imaging and quantification 

Automated imaging was carried out using the Cell Lab 
IC-100 Image Cytometer (IC100, Beckman Coulter) at 
room temperature using CytoShop 2.0. Image acquisition 
was performed with a Nikon S Fluor 40 x /0.90 NA object- 
ive. Z-stacks were imaged at 1 uM intervals at 1 x 1 
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binning. Nuclear array segmentation and signal quantiza- 
tion were performed using PipelinePilot image analysis 
software as previously described (20). Briefly, maximum 
intensity projections were created for the GFP-ERa and 
antibody channels followed by background subtraction to 
all images. Nuclei were defined by DAPI signal and masks 
created by applying adaptive thresholding followed by 
marker-based watershed. Array segmentation was based 
on thresholding of the GFP-ERa image. Aggregated cells, 
mitotic cells and apoptotic cells were removed using filters 
based on nuclear size, nuclear shape and nuclear intensity. 

Deconvolved images were taken using a DeltaVision 
Core Image Restoration Microscope (Applied Precision) 
at room temperature. Z-stacks were imaged at 0.2 uM 
intervals with a 60x/1.42NA Plan Achromat objective 
and a charge couple device camera (CoolSnap HQ2; 
Photometries). These images were also taken at 1 x 1 
binning. 

Quantitative reverse transcriptase-polymerase chain 
reaction and chromatin immunoprecipitation 

MCF-7 cells were treated with 10 nM E2, 100 nM Dex or 
E2 + Dex for 45 min before the chromatin immunopre- 
cipitation (ChlP) procedures. Cells were fixed with 1% 
formaldehyde for 10 min at room temperature, lysed in 
RIPA (Sigma-Aldrich) buffer, sonicated and then the 
soluble fraction pre-cleared for 1 h with protein A/G 
agarose beads. Pull-downs were then performed using 
the HC-20 anti-ERa antibody (Santa Cruz) and Pierce 
Rbt anti-GR (Thermo Scientific). After the IP, protein 
A/G agarose beads were added and washed sequentially 
with high salt, low salt and TE buffers. DNA was then 
eluted and analyzed by quantitative real-time polymerase 
chain reaction (PCR) using specific primers for the 
ER-binding sites closest to selected ER-regulated genes. 

Total RNA was isolated using Trizol® as per manufac- 
turer's instructions. One microgram of total RNA was 
reverse transcribed in a total volume of 20 ul using 200 
U reverse transcriptase, 50pmol random hexamers and 
1 mM dNTPs (New England Biolabs, Beverly, MA, 
USA). The resulting cDNA was analyzed with real-time 
PCR. Each real-time PCR reaction consisted of 6ul of 
diluted reverse transcribed product, lx SYBR Green 
PCR Master Mix (Applied Biosystems, Foster City, CA, 
USA) and 50 nM of forward and reverse primers. 
Reactions were carried out in an Applied Biosystems 
StepOne Plus for 40 cycles (95°C for 15 s, 60°C for 
1 min). The fold change in expression was calculated 
using the A A Ct comparative threshold cycle method, 
with the ribosomal protein 36B4 mRNA as an internal 
control. 

Statistical analysis and software 

Data presented were acquired from a minimum of three 
independent experiments performed on different days with 
different passages of cells. P-value was determined using 
student's ?-test for single comparisons or one-way 
ANOVA for independent samples and a Tukey HSD 
post-hoc test for comparisons across multiple samples. 



The data for the heat map in Figure 5A were generated 
using Cluster 3.0 using hierarchical clustering with a 
centroid linkage method (28) and then visualized using 
Java Treeview (29). 

RESULTS 

Type I nuclear receptors load onto the PRL array only 
after ERa recruitment 

We have previously generated and characterized the 
PRL-HeLa array cell line stably expressing GFP-ERa 
(GFP-ERa:PRL-HeLa) as a model to visually study and 
inspect the mechanisms linked to ERa mediated transcrip- 
tional regulation by high-throughput fluorescence micros- 
copy and high-content analysis (19-22,30,31). To expand 
our knowledge of the interplay between multiple tran- 
scription factors at the integrated transcription locus, we 
performed bioinformatic analysis of the promoter/ 
enhancer elements that comprise the PRL array to deter- 
mine the full complexity of additional transcription factor- 
binding sites. Using web-based motif analysis software 
(e.g. Alibaba and JASPAR), we identified several 
putative transcription factor- binding motifs (CEBP, API 
and GATA) in the PRL promoter/enhancer, including a 
steroid/hormone response element (HRE/SRE) located 
close to the ERE/Pitl synergy element (Figure 1A for 
diagram). 

With our longstanding interest in the molecular biology 
of type I NRs coupled with the recent findings from the 
various laboratories on transcription factor interplay (12- 
15), we asked whether there were endocrine conditions 
where the HREs/SREs within the PRL array could be 
occupied by combinations of type I NRs. Although all 
GFP-ERa:PRL-HeLa cells exhibit nuclear ERa, the 
multicopy PRL locus becomes visibly 'loaded' only when 
ligand treated; when it is grown in media containing 
stripped/dialyzed serum, virtually no nuclei exhibit a 
visible array. To this end, we transfected the parental 
PRL array-containing HeLa cells (PRL-HeLa, which are 
ERa negative) with GFP-tagged versions of five type I 
NRs e.g. AR, PR-A, PR-B and GR) either alone (not 
shown) or in co-transfection with mCherry-ERa and 
treated them for 30 min with lOOnM dexamethasone 
(for GR), 10 nM dihydroxytestosterone (DHT) (for AR) 
or 10 nM R5020 (for PR-A and PR-B). Our results are 
summarized in Figure 1B-E and reveal that under NR 
agonist alone conditions (except for E2, see later in the 
text) in stripped/dialyzed serum-containing media, 
GFP-AR (Figure IB), GFP-GR (Figure 1C), GFP-PRB 
(Figure ID) and GFP-PRA (Supplementary Figure SIB) 
all failed to load onto the PRL array. Further, as 
expected, the non-ER type I NRs did not load the PRL 
array with E2 treatment; in contrast, 17P-estradiol (E2) 
treatment readily induced array formation in all ERa-ex- 
pressing cells (middle column Figure 1B-E). However, 
combining E2 treatment with any of the agonists for 
type I NRs (i.e. dexamethasone, DHT and R5020) 
induced loading of both ERa and type I NRs as high- 
lighted by the arrows (third column Figure 1B-E) and is 
indicative of an ERa 'cooperative loading' mechanism. As 
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Figure 1. Type I NRs localize to the PRL array only under dual ligand conditions. (A) Schematic representation of the PRL array that was stably 
integrated into the PRL-HeLa array cell line along with GFP-ERot. (B) Images of transfected ERa (bottom row) and AR (top row) in PRL HeLa 
cells treated with lOnM DHT, lOnM E2 or E2 + DHT. (C) Images of transfected ERot (bottom row) and GR (top row) in PRL HeLa cells treated 
with 100 nM Dex, lOnM E2 or E2 + Dex. (D) Images of transfected ERa (bottom row) and PR-B (top row) in PRL HeLa cells treated with lOnM 
R5020, lOnM E2 or E2 + R5020. (E). Images of transfected ERoi (bottom row) and PPARy (top row) in PRL-HeLa cells treated with lOnM Rosi, 
lOnM E2 or E2 + Rosi. White arrows indicate the location of NR loading onto the PRL array. 



an additional control, we also examined a non-type I NR, 
PPARy (Figure IE), which does not recognize HRE/SRE- 
binding motifs. Indeed, PPARy was unable to load onto 
the PRL array either after treatment with its own agonist 
(rosiglitazone) alone or in combination with E2 (Figure 
IE), suggesting that the 'cooperative loading 1 observed in 
Figure 1 does not occur with all NRs. 

We were also interested in whether ERp, which can 
localize to the PRL array with E2 (31), was able to 
cause recruitment of the type 1 NRs to the PRL array. 
To determine whether the cooperative loading was able to 
occur with ERp, we transfected GFP-GR and ER|3- 
mCherry and treated with lOnM E2+100nM Dex. As 
seen in Supplementary Figure SIB, GR was able to 
localize to the PRL array element with ERp. Because of 
the relatively low homology between ERa and ERp 
N-termini, these data are also suggestive that the 
N-terminus of ERa may not be important for the coopera- 
tive loading. 

To determine whether the NRs binding to the PRL 
array is dependent on DNA binding to the SRE/HRE, 
we transfected wild-type PR-B or PR-B with a mutated 
DNA binding domain (DBD), which makes it unable to 
bind directly to DNA, PR-B DBDmut (32), into 
PRL-HeLa cells constitutively expressing GFP-ERa 
and treated them with Veh or 10 nM E2+ lOnM R5020. 



Wild-type PR-B was able to colocalize to the PRL array 
with ERa on dual agonist treatment (Supplementary 
Figure SIC), whereas the PR-B DBDmut was unable to 
do so under either condition (Supplementary Figure SID). 
These data strongly suggest that localization to the PRL 
array of the type 1 NRs after treatment with E2 + agonist 
is dependent on DNA binding to a SRE/HRE sequence. 

Loading of endogenous GR onto the PRL array is 
ERa-dependent and maximal when both receptors 
are agonist bound 

Next, we wanted to confirm that the observed 'cooperative 
loading' between NRs was not due to overexpression of 
the partners or the presence of tags (GFP and mCherry). 
We decided to focus on GR due to its endogenous expres- 
sion in the PRL-HeLa array cell line, and due to the 
growing body of literature discussing different modes of 
cross-talk between glucocorticoids and estrogens (33-35). 
To assess loading of endogenous GR onto the PRL array, 
we performed antibody labeling for GR in the 
GFP-ERa:PRL-HeLa cell line treated with vehicle, 
10 nM E2, 100 nM Dex or E2 + Dex for 30 min. As 
expected, GR translocated into the nucleus (Figure 2A, 
bottom panels) under Dex or E2 + Dex treatments, but 
GR loading onto the PRL array was only evident after 
dual agonist treatment (E2 + Dex; arrows in Figure 2A as 
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Figure 2. Endogenous GR localizes to the PRL array under dual agonist conditions. (A) Representative fields of PRL-HeLa:GFP-ERa cells treated 
with 5% EtOH (veh), 10 nM E2, 100 nM Dex or 10 nM E2+ lOOnM Dex, showing DAPI, GFP-ERot and GR staining. (B) Quantification of ERa 
and GR loading index (mean intensity of signal at array/mean intensity of signal in the nucleoplasm) using PipelinePilot Image Analysis software. 
*P<0.05. More than 1000 cells/condition were analyzed. (C) ERa and GR were transfected into PRL-HeLa cells treated with E2 + Dex and show 
localization to the PRL array. (D) mCherry-tagged ERa and GFP-tagged GR were transfected into PRL-HeLa cells and treated with E2 + Dex. Only 
the cell that has both ERa and GR shows GR localization to the PRL array (white arrows). 



examples). This was also true when ERp, instead of ERa, - 
containing cells were used (Supplementary Figure SIB). 
We quantified the amount of GR that loaded onto the 
PRL array (Figure 2B) by using a custom image analysis 
protocol developed on the Pipeline Pilot Image Analysis 
platform (PLP; see 'Materials and Methods' section for 
details). The PLP platform permits us to perform robust 
analysis of ERa functions by segmenting array and nu- 
cleoplasm regions and then quantify, on a cell-by-cell 
basis, the amount of any labeled factor present at the 
PRL array. We calculated a loading index for GR and 
ERa, which is defined as the mean intensity of GR or 
ERa signal at the array divided by the mean intensity of 
GR or ERa signal in the nucleoplasm. Using this metric 
with automated microscopy, we were able to quantify 
> 1000 cells per treatment and observed a significant 
increase in ERa loading to the PRL array on E2 and 
E2 + Dex treatments (with no effect of Dex on ERa 
loading); in contrast, GR loading onto the PRL array 



was significantly increased only after E2 + Dex treatment 
(Figure 2B). We also performed time course and hormone 
dose-response curves and measured ERa and GR loading, 
and per cent of cells with an array and array size to fully 
characterize this cooperative interplay (Supplementary 
Figure S2). 

We then wanted to determine whether the loading of 
GR onto the PRL array was influenced by the GFP-tag 
fused to ERa. To verify this was not the case, we trans- 
fected untagged ERa or mCherry-ERa into parental 
PRL-HeLa cells and treated them with E2 + Dex to deter- 
mine whether GR still exhibited cooperative loading. 
After dual agonist treatments, GR localized onto the 
array with each version of ERa demonstrating that the 
cooperative loading is not dependent on the type, 
presence or location of the fluorescent tag. The 
immunolabeling images for untagged ERa are shown in 
Figure 2C. To discern whether GR loading onto the PRL 
array is ERa-dependent or may be due to other signaling 
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pathways (e.g. via non-genomic actions or GPR30 activa- 
tion), we co-transfected PRL-HeLa with GFP-GR and 
mCherry-ERoc and treated them with E2 + Dex for 30 
min. Figure 2D shows that GR recruited onto the PRL 
array only in dual-transfected cells (white arrow). In cells 
where no ERoc expression is evident, GR did not target the 
PRL array, indicating that cooperative loading is depend- 
ent on nuclear ERa loading onto the PRL array. 

We next wanted to evaluate whether cooperative loading 
was dependent on specific classes of ligands (i.e. agonists 
versus antagonists) for either ERa or GR. To this goal, 
we first treated the GFP-ERa:PRL-HeLa cells with 
vehicle, 100 nM Dex or the GR antagonist RU486 
(lOOnM) in the presence of 10 nM E2 for 30 min (Figure 
3A). Although RU486 treatment did cause an increased 
GR loading onto the PRL array, it was significantly less 



compared with Dex treatment (33% reduction). Similarly, 
when we treated with ERa antagonists [10 nM 
4-hydroxy-tamoxifen (Tam) or 100 nM Fulvestrant (ICI)] 
in the presence of lOOnM dexamethasone (Figure 3B), we 
saw that both Tam and ICI were able to cause recruitment 
of GR onto the PRL array but at a significantly lower level 
than E2 + Dex (66% reduction). Representative images are 
shown in Supplementary Figure S3. Collectively, the data 
shown in Figure 3A and B suggest that GR loading onto 
the PRL array is maximal when both receptors are in an 
agonist conformation. 

GR loading onto the PRL array correlates with ERa 
loading but not array area, GR level or ERa level 

High content analysis allows us to mine image data sets 
with hundreds of features per individual cell. These 
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Figure 3. ERa, E2 and Dex are required for maximal GR loading. (A) Quantification of GR loading index in the presence of 10 nM E2 with either 
5% EtOH (veh), 100 nM Dex or 100 nM RU486. (B) Quantification of GR loading index in the presence of 100 nM Dex with vehicle, 10 nM E2, 
10 nM Tam or 100 nM ICI. *P<0.05 compared with veh treatment. (C-F) Single-cell scatter plot of GR loading index (.\-axis) with ERol loading 
index (C), array area (D), ERa nuclear intensity (E) or GR nuclear intensity (F) on the v-axis. Values above each plot are correlation coefficients 
quantified using Pearson's product moment. 
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features are descriptors of various mechanistic or cytolo- 
gical characteristics of the fluorescent signals being 
examined; they are linked to protein expression, intracel- 
lular and intranuclear localization and so forth. We chose 
four basic features directly linked to transcriptional regu- 
lation, ERa binding to DNA, chromatin remodeling, ERa 
nuclear protein level and GR nuclear protein level, and we 
used them to determine whether a correlation exists 
between this feature set and GR loading onto the PRL 
array. Figure 3C-F show scatter plots of single cell 
values for E2 + Dex treatment with the x-axis being the 
GR loading index and j-axis corresponding to one of the 
four chosen features. No statistically significant correl- 
ation was observed between GR loading index and array 
area (Figure 3D), ERa nuclear levels (Figure 3E) or GR 
nuclear levels (Figure 3F). Conversely, GR loading index 
did significantly correlate (Pearson's r = 0.5126) with ERa 
loading index (Figure 3C), providing additional evidence 
that an ERa loaded PRL array is required for GR PRL 
array occupancy. 

GR loading onto the PRL array is dependent on 
ERa helix-12 

To further explore the mechanism(s) through which ERa 
assists GR recruitment onto the PRL array, we wanted to 
determine which domain of ERa is required. We transi- 
ently co-transfected PRL-HeLa cells with mCherry-GR 
and GFP-tagged wild-type ERa or mutant ERa vectors 
that have previously been used in the study for the regu- 
lation of ERa (36). We compared wild-type ERa with a 
helix-12 truncation (amino acids 1-534), which eliminates 
the interaction domain for the pi 60 family of co-regula- 
tors (e.g. SRC-3) and an F domain truncation (amino 
acids 1-554) (Figure 4A). Cells were then treated with 
E2 + Dex for 30 min and immunolabeled with an 
antibody to RNA polymerase II (Pol II) and GR 
(Figure 4A), or SRC-3 and GR (Figure 4B). As 



expected, wild-type ERa was able to recruit GR, Pol II 
and SRC-3. The F-domain mutant (amino acids 1-554) 
was also able to recruit GR, Pol II and SRC-3. The 
mutant ERa lacking helix-12 (amino acids 1-534) was 
capable of loading onto the PRL array; however, it 
failed to recruit GR, Pol II or SRC-3. These data 
indicate that GR recruitment onto the PRL array is 
mediated through helix-12 and requires an intact ligand- 
binding domain and, as shown later in the text in detail, 
may involve other co-regulators (Figure 5). 

MED14, SRC-2 and SRC-3 mediate GR loading onto 
the PRL array 

To identify potential co-regulators that mediate GR 
loading onto the PRL array, we performed a small 
siRNA survey directed against a panel of well-established 
modulators of transcription, including the pi 60 family of 
co-regulators (SRC-1, -2 and -3), LEOl, DDX3X, 
MED 14 and Sin3B. We also used GR siRNA as positive 
control and scrambled siRNA as negative control. 
GFP-ERa:PRL-HeLa cells were transfected with the 
siRNA panel for 72 h and were then treated for 30 min 
with 10 nM E2+100nM Dex; after fixation, cells were 
immunolabeled to localize GR and the more active 
serine 5-phosphorylated form of RNA polymerase II. 
We were thus able to perform high-content analysis to 
determine the effects of each siRNA on nuclear intensity 
(levels), array intensity (array int) and array loading for 
ERa, GR and serine 5 phosphorylated Pol II. Figure 5A 
shows a heat map summarizing the results from the 
siRNA knockdown analysis. Although most of the genes 
tested (SRC-1, SRC-2, SRC-3 and so forth) have little 
effect on PRL array features (SRC-3 siRNA results in a 
loss of active Pol II at the array, whereas Sin3B has the 
opposite effect), the three outliers (combination SRC-2/-3, 
GR and MED 14) all cause a significant decrease in GR 
loading. Light blue indicates reduction, whereas yellow 
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Figure 4. ERa mutants demonstrate importance of ERa helix-12 in GR recruitment. (A and B) GFP-tagged ERa mutants (full length, amino acids 
1-534 or amino acids 1-554) were transfected along with mCherry-tagged GR into PRL-HeLa cells and treated with E2 + Dex for 30 min. Cells were 
labeled with antibody targeting RNA polymerase II (A) or SRC-3 (B). White arrows indicate the location of NR, SRC-3 or RNA polymerase II 
loading onto the PRL array. 
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Figure 5. MED 14, SRC-2 and SRC-3 mediate the recruitment of GR to the PRL array. (A) Heat map showing effects of siRNA knockdown of 
multiple factors on features of the PRL array system, including GR loading, active RNA polymerase II loading and ER loading. (B) Quantification 
of GR nuclear intensity (normalized to control siRNA) after 72h siRNA knockdown of GR, MED14 and SRC-2/-3 in HeLa-PRL:GFP-ERa cells 
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after 72h siRNA knockdown of MED 14, or SRC-2/-3 in HeLa-PRL:GFP-ERa cells treated with lOnM E2+100nM Dex for 30min (*,P<0.05; 
**p< 0.005; *** P< 0.0005). 



reflects an increase in each of the nine parameters 
examined compared with scrambled siRNA (black). As a 
control, in Figure 5B, we show effective GR knockdown 
(e.g. reduced GR nuclear intensity and array loading 
without affecting ER levels or array loading). 

We next tested which co-regulator knockdowns 
specifically impacted GR loading on the PRL array 
without affecting global GR nuclear levels. High-content 
analysis of the RNAi survey identified MED 14, a subunit 
of the mediator complex, and to a lesser extent SRC-2, 



as essential for GR PRL array loading after treatment 
with both E2 and Dex (Figure 5A). Intriguingly, a com- 
bination of SRC-2 and SRC-3 siRNAs produced a greater 
reduction in GR loading (Figure 5E) compared with single 
knockdowns of the pl60 co-regulators. We validated 
siRNA knockdown efficacy by immunofluores- 
cence microscopy as shown in Figure 5B-D and 
Supplementary Figures S4 and S5D. SRC-2, SRC-3 
and MED 14 are all recruited to the PRL array on both 
E2 and E2 + Dex conditions with no statistical difference 
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between the two treatments (Supplementary Figure 
S5A-C). 

To infer the spatial interrelationships between SRCs 
and MED14 that lead to GR loading, we knocked down 
GR (Figure 5B), MED14 (Figure 5C) and SRC-2/-3 
(Figure 5D-E, respectively) and then quantified the 
loading of GR (Figure 5F), MED 14 (Figure 5G), SRC-2 
(Figure 5H) and SRC-3 (Figure 51). Knockdown of 
SRC-2/-3 decreased loading of GR and MED14, and, as 
expected, also SRC-2, and SRC-3, whereas MED 14 
siRNA only blocked MED14 and GR loading. These 
data suggest that SRC-2 and SRC-3 are required for 
MED 14 loading, which is then responsible for GR 
loading. 

GR presence at the PRL array causes reduction in 
E2-induced reporter gene transcription 

To better understand the functional consequences of ERa- 
GR cross-talk, we sought to directly quantify the effect of 
E2 + Dex treatment directly on transcriptional output 
from the PRL array by using RNA FISH directed 
against the dsRED2 reporter gene mRNA. GFP- 
ERa:PRL-HeLa cells were treated with vehicle, E2 or 
E2 + Dex for a time course ranging from 5min to 24 h, 
and accumulation of the PRL-regulated dsRED2 
reporter gene mRNA was measured by RNA FISH 
using a highly specific fluorescently labeled probe set 
(Stellaris T , Biosearch Technologies). As shown in 
Figure 6A, the dsRED2 mRNA was rapidly induced by 
E2 treatment and E2 + Dex treatment alike. Around 6h 
after treatment, E2 + Dex treatment resulted in signifi- 
cantly less mRNA than E2 treatment alone, indicating 
that loading of GR on the PRL array ultimately leads to 
a transcriptional repressive effect. 

Complex ERa/GR interplay on chromatin in MCF-7 
breast cancer cells 

We next wanted to determine whether ERa/GR coopera- 
tive chromatin loading and transcriptional interplay was 
occurring at any known ERa target genes in an endogen- 
ous setting. Using MCF-7 breast cancer cells that had 
been treated with 5% EtOH, 10 nM E2, 100 nM Dex or 
E2 + Dex for 24 h, we performed qPCR on two well- 
characterized ERa target genes, TFF1 (pS2, Figure 6B) 
and GREB1 (Figure 6D). TFF1 mRNA levels, increased 
with E2, were reduced with Dex and showed a muted E2 
response in the E2 + Dex treatment. We then performed 
ChIP assays for ERa and GR to determine their recruit- 
ment to the TFF1 enhancer after 45 min of treatment 
(Figure 6C). As expected, ERa recruitment to the TFF1 
enhancer was increased by E2 treatment and did not 
change by addition of Dex. We also show that GR is re- 
cruited to the same enhancer under both Dex and 
E2 + Dex treatments without showing any statistically sig- 
nificant difference. In contrast to TFF1, GREB1 is 
induced by both E2 and Dex with a synergistic increase 
in the double treatment. On analysis of ERa ChlP-seq 
data (37), it is clear that the GREB1 genomic locus is 
rich in ERa binding sites. By ChlP-PCR analysis, we 
found two of these ER-binding sites (GREB_732 and 



GREB_734) that show significant increases in GR 
binding compared with vehicle (Figure 6E and F, respect- 
ively), but only when treated with E2 + Dex, with minimal 
effects on ERa binding (Figure 6E-F). The effect of Dex 
alone is likely because of novel GR-binding elements in 
the GREB1 locus, as indicated by ChlP-seq analysis (K. P. 
White, personal communication). These data not only re- 
capitulate what we observed in the PRL array model 
system but also highlight gene- and enhancer-specific 
effects that can occur in an endogenous setting, which 
we are continuing to explore. 



DISCUSSION 

NRs are a key group of transcription factors that regulate 
gene regulation programs in many diverse tissues where 
they serve as rheostats that respond to external and 
hormonal stimuli. The organismic distribution of NRs 
shows both widespread and tissue-specific expression 
(38), indicating that extensive interplay between 
numerous hormones at the same time must occur. 
Moreover, it is well known that combinatorial expression 
of NRs in tumors (i.e. breast cancer) is important in 
determining prognosis and avenues of treatment (39,40). 
Despite the obvious complexity of multiple NRs function- 
ing simultaneously in any cellular milieu, thus far studies 
interrogating NR function(s) have been most extensively 
examining only one hormone at a time; regardless, these 
highly reductionist approaches have provided early essen- 
tial insight into NR actions. Notable exceptions include 
studies comparing GR with PR actions (41,42) or ERa 
with retinoic acid receptors (12,13); collectively, these 
efforts showed both competitive and cooperative inter- 
actions. To quantitatively explore the functional complex- 
ity and interplay of type I NRs with visual tools, we have 
used an engineered multicopy estrogen responsive tran- 
scriptional array system (PRL array) and high-throughput 
microscopy as an assay platform. The PRL array system 
affords us with a unique opportunity to visualize by fluor- 
escence microscopy and quantify, on a cell-by-cell basis 
using high-content analysis, all the major steps involved 
in ERa-mediated transcriptional activation, including 
receptor and co-regulator recruitment, histone modifica- 
tions, chromatin changes and transcriptional output. 

After bioinformatics analysis that revealed an add- 
itional NR-binding site in the PRL enhancer, we used 
the PRL array model to explore potential functional inter- 
play between ER and other type I NRs at the chromatin 
level. Our observations that dual ligand treatments (e.g. 
E2 + another agonist) are required for ERa to recruit GR 
(or other type I NRs, but not PPARy) to the adjacent 
HRE/SRE are consistent with an assisted loading mech- 
anism similar to recently described by Voss et al. (15). 
However, Voss et al. analyzed the cross-talk of two NRs 
on a single ra-element using wild-type GR and an ERa 
point mutant that alters ERa specificity from ERE to 
GRE binding. In the current study, we used the PRL 
model to analyze the interplay between two distinct cis- 
elements (ERE versus HRE/GRE) using two wild-type 
NRs (ERa and GR). 
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Figure 6. ERa/GR cross-talk effects transcriptional output in PRL array and MCF-7 cell systems. (A) Time course of dsRED2 mRNA production 
at the PRL array when treated with 10 nM E2 or 10 nM E2+ lOOnM Dex and measured by RNA FISH. *,P<0.05 versus E2. (B) TFF1 mRNA 
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normalized to vehicle treatment. (D) GREB1 mRNA levels were measured by qPCR after 24h of treatment. (E and F) ChIP assay showing ERa 
(gray) and GR (black) recruitment to two enhancer regions (732 and 734, respectively) of the GREB1 genomic locus after 45-min treatments. 



We report here our quantitative observations of de novo 
loading of type I NRs onto the PRL array only under dual 
ligand treatment (i.e. ERa ligand + type I NR ligand). 
Collectively, these data suggest a central role for ERa-de- 
pendent modulation of high-level chromatin structure that 
permits other transcription factors to access their respect- 
ive binding elements that were previously unavailable. We 
were also able to show, using a PR-B DBD mutant, that 
this localization requires DNA binding. Additionally, we 



used this model system to illustrate the cooperative 
loading phenomenon is specific for type I NRs, as 
PPARy did not localize to the PRL array under dual 
ligand treatment. We cannot, however, exclude that 
members of other families of transcription factors could 
follow the same mechanism in this model, or that other 
NRs could use the same mechanism on a selected number 
of response elements in the whole genome. Studies on 
these scenarios are currently under way, but combined 
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with the work of Voss et al. is suggestive that type I NRs 
have this pioneering ability. 

Our decision to further explore the mechanisms 
underlying the interplay between ERa and GR at the 
chromatin level is emphasized by the importance of both 
receptors in prognosis and treatment of breast cancer (40) 
and other pathologies including osteoporosis and inflam- 
matory diseases (43). Moreover, GR is the only type I NR 
that is endogenously expressed in HeLa cells, thus 
facilitating mechanistic studies without the need of ex- 
ogenous expression. Our development and application of 
custom, single-cell-based high-content analysis allowed us 
to quantify multiple mechanistic steps that are involved 
during or after de novo loading of endogenous GR to 
the PRL array. To address which domain of ERa was 
essential for recruitment of GR under dual ligand treat- 
ment, we used a series of well-characterized ERa mutants, 
and we found that helix- 12 within the ligand-binding 
domain of ERa is required for GR loading onto the 
PRL array. Although other ERa domains may be 
playing a role as well (e.g. the DNA-binding domain), 
we observed a complete loss of GR loading with the 
helix- 12 mutant, indicating a central role for this portion 
of the receptor. Moreover, ER(3, which has a different 
N-terminal domain from ERa, is still fully capable of re- 
cruiting GR. As helix- 12 is the main structural feature that 
determines recruitment of co-regulatory factors, this ob- 
servation raised the possibility that coactivators and their 
chromatin remodeling actions may be necessary require- 
ment for GR recruitment to the array. 

To test this hypothesis, we performed a small focused 
siRNA survey against selected transcriptional modulators. 
Our survey identified MED 14, SRC-2 and SRC-3 as three 
co-regulators that are required for GR loading after dual 
agonist treatment, without effecting changes on ERa 
loading to the PRL array. Interestingly, a reduction in 
GR loading also occurred when we treated PRL- 
HeLa-ERa cells with E2 + RU486 after knockdown of 
MED14, suggesting a similar method of recruitment re- 
gardless of GR ligand (data not shown). This was an inter- 
esting observation, as MED14 and both SRC-2 and 
SRC-3 have been previously studied in the context of 
GR actions. MED14 has been shown to interact with 
the N-terminal AF-1 region of GR (44), leaving the 
C-terminal AF-2 region free to interact with other tran- 
scriptional regulators. In breast cancer cell lines, MED14 
has also been shown to act as an ERa co-activator 
through helix- 12 interaction (45), further implicating this 
component of the mediator complex as a scaffold for re- 
cruitment of additional factors. This is the first data to 
suggest a novel integrator role for mediator proteins to 
modulate recruitment to sites throughout the genome of 
different transcription factors under various hormonal 
stimuli. 

From our initial RNAi survey, we also found that 
SRC-2 and SRC-3 are required for MED 14 recruitment 
to the PRL array, but not vice versa, indicating that 
MED14 may be recruited as a second tier co-regulator 
to help recruit GR. Our study is also the first data on 
coactivators playing an important role in the assisted 
loading of transcription factors. It is also interesting to 



observe the compensatory effects of SRC-2 and SRC-3 
on the ability of GR to load onto the PRL array. 
Previous reports have shown fewer compensatory (46,47) 
than non-compensatory (48-50) interactions between 
pl60 family of coactivators. These suggest a greater role 
in complex transcriptional mechanisms for the pi 60 
family members beyond transactivation and suggest modi- 
fication of cross-talk in cancers that overexpress the 
coactivators. 

The effect of the ERa/GR interplay was also queried at 
the level of transcriptional regulation, and we found that 
the presence of both GR and ERa at the PRL array causes 
a reduction in E2-mediated stimulation of gene transcrip- 
tion. However, this was not the only possible outcome, as, 
in MCF-7 cells, we observed both cooperation and antag- 
onism on E2 regulated genes, as exemplified by the TFF1 
and GREB1 mRNAs. The mechanism(s) underlying this 
phenomenon is not fully clear, and it is likely to be gene 
specific, with both issues calling for future investigations. 
We were also able to show cooperative interplay between 
ERa and GR in MCF-7 cells with a distinct transcrip- 
tional output at the GREB1 and TFF1 loci. The different 
styles of interplay in the enhancer regions of GREB1 and 
TFF1 demonstrate the need for additional studies 
examining how transcription factors are arrayed along 
the genome when multiple hormones are present, as 
clearly this is the more physiologically relevant context. 
Taken together, the results presented here highlight a 
mechanism by which GR (and/or other steroid NR recep- 
tors) can be recruited to previously inaccessible chromatin 
regions by a 'pioneering' role of ERa in chromatin re- 
modeling and co-activator interaction that ultimately 
modulates ERa-regulated gene transcription. In particu- 
lar, GR requires the combination of ERa, MED14 and 
SRC-2/-3 to be recruited. However, we and others have 
failed to observe a stable complex formation between ERa 
and GR using both co-immunoprecipitation and ChlP- 
reChlP experiments [data not shown and (15)], suggesting 
a dynamic relationship between transcription factors and 
co-regulators at enhancer elements. The observation of 
this interplay is an important step towards further dissec- 
tion of steroid-mediated cross-talk between pathways that 
underlie transcriptional fine tuning to meet cellular re- 
quirements for homeostasis, or play aberrant roles 
during pathogenesis and disease development. 
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